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A time-dependent Synchrotron Self Compton model (SSC) which is able to motivate 
the used electron spectra of many SSC models as a balance of acceleration and radia- 
^ ' tive losses is introduced. Using stochastic acceleration as well as Fermi-I processes even 

(J"^ electron spectra with a rising part can be explained, which are mandatory to fit the 

^y\ , lowstate spectral energy distribution (SED) of PKS 2155-304 as constrained from Fermi 

^^ ' LAT observations. Due to the time resolution the outburst of PKS 2155-304 observed 

^^ ' by H.E.S.S. in 2006 can be modelled selfconsistently as fluctuations along the jet axis 

without introducing new sets of parameters. The model makes the time evolution of the 
SED also accessible. Hence giving new insights into the flaring behavior of blazars. 



Keywords: jets; radiation mechanisms: non-thermal; BL Lacertae objects: PKS 2155-304 

1. Introduction 



Blazars have gained a lot of attraction over the last years mainly due to new obser- 
j^ ■ vations from imaging air Cherenkov telescopes like H.E.S.S., MAGIC or VERITAS 

and lately through the Fermi satellite. It is now possible to cover major parts of the 
very high energy (VHE) band from 30 GeV up to 30 TeV. The nonthermal emission 
of blazars is assumed to be Doppler enhanced radiation from the relativistic jet of 
an Active Galactic Nuclei (AGN) which axis is covering only a small angle with the 
line of sight. The spectral energy distribution (SED) of blazars exhibits two distinct 
humps. Depending on the position of those maxima blazars are divided into at least 
two subcategories, the highly luminous Flat Spectrum Radio Quasars (FSRQs) with 
the first peak occurring in the optical regime and the BL Lac Objects (high and 
low peaked) having the first hump up to the hard X-Rays. 

The cause of the nonthermal radiation is not yet clarified, but a number of models 
have been developed and successfully applied to single objects such as hadronic 
models Jl|, Synchrotron Self Compton (SSC) and External Compton (EC) models 
[e.g. I2, y or the lately adopted spine layer ones [e.g. l4i]. Unlike FSRQs BL Lac 
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objects seem to be dominated by leptons, therefore the Compton models, where the 
first hump in the SED is due to synchrotron radiation and the second to upscatter- 
ing of photons by electrons are very successful. We will concentrate on a SSC model 
where internal photons provide the targets for the inverse Compton effect. Even 
though these models are commonly used, there are still open questions. The first 
one our model will try to answer concerns the electron distribution, most models 
assume a fixed powerlaw rather than applying a selfconsistent evolution of electron 
spectra. The second one is, how the variability may be explained in this context. 
Standard models are not time dependent and may therefore not model variability. 
The model we developed is a timede-pendent SSC model, which explains the elec- 
tron distribution as a balance of acceleration and radiation. 

To test our model we choose PKS 2155-304 (henceforth shortly PKS2155) a HBL 
at redshift z = 0.117 i) for it is well observed over a long period of time with the 
H.E.S.S. telescope giving a high confidence on the quiescent state, ii) for the re- 
cent multiwavelength campaign including Fermi l5l| giving a good coverage of the 
whole SED of PKS2155. iii) PKS2155 is a very variable source with its famous 
fiare in 2006 [6|, the strongest outburst ever measured at these VHE. All these 
circumstances making it the most suitable blazar for our investigations. 

2. The Model 

In our model we solve the time-dependent kinetic equation, derived from the one 
dimensional diffusion approximation [e.g.^] of the Vlasov equation [see e.g.iSj using 
the relativistic approximation p w jmc, in two spatially different zones. Both zones, 
the acceleration zone and the radiation zone, are assumed to be spherical and ho- 
mogeneous containing isotropically distributed electrons and a randomly oriented 
magnetic field. All calculations are made in the rest frame of the blob. Hence elec- 
trons enter the spatially smaller acceleration zone (radius i?acc) upstream the blob 
and are continuously accelerated due to shock and stochastic acceleration being bal- 
anced by synchrotron radiation, this extends the model of [9|]- After tosc = V^ncc/c 
the electrons enter the radiation zone (radius i?rad) where rj is an empirical fac- 
tor set to 7] = 10. Downstream the electrons are suffering synchrotron as well as 
inverse Compton losses. Other contributions like pair production do not alter the 
SED in typical SSC conditions and are neglected [3[. To calculate the model SED 
these photons arc Doppler shifted towards the observer taking the redshift z of the 
source into account, i.e. 1,^^^^^ = i5'^/ifobs/(47r)A^ph with fobs = <5/(l -I- z)u. Due to 
its size the acceleration zone does not contribute to 7obs directly, since the observed 
flux at distance r is given as -Fi/„b, (r) — Ti'-^i^oba-^rad'^^^- ^he kinetic equation in the 
acceleration zone is 



—Q^ =^ mi - t,..l) ■ neh, t)\ + — 



[(a-f 2)tacc]"'7 



-1 2dneh,t) 



dj 



+ Qo{l-lo)-tcslneil,t) . (1) 
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The injected electrons, when the blob propagates through the jet, are considered 
via Qinj(7,i) := QoS{j — 7o)- The synchrotron losses are calculated using eq. ([2]). 

Psh) - ^^^l' - M' (2) 

OTT mc 

with the Thomson cross section ctt- One finds iacc — («s/(4^||) + '^v\/{9K\\)^ 
for the characteristic acceleration tinicscale of the system by comparing eq. ([1]) 
with [7|. K\\ is the parallel spatial diffusion coefficient which does not depend on 
7 when using hard spheres. Furthermore, unlike in [10| . the energy dependence of 
the escape losses is also neglected since we do not expect a pileup as suggested in 
[7| at typical SSC conditions. Vs.va are the shock and Alfven speed respectively, 
hence a in eq. ((!]) is proportional to v1/v\. All escaping electrons enter the radiation 
zone downstream the jet. This ansatz takes account of a much more confined shock 
region, for Fermi-I acceleration will probably not occur over the whole blob when 
considering real blazars. 
The lack of acceleration in the radiation zone simplifies the kinetic equation to 

OT 07 trad,osc \ i?rad / tosc 

The additional Pic accounts for the inverse Compton losses of the electrons occur- 
ring beside the synchrotron losses [e.g. |8|: 

Pic(7) = m^c'h r^^^ daa H da,N,,,,{a,) '^^^]^'^'^ (4) 

Jo ^0 dtda 

With the photon energies rewritten in terms of the electron rest mass, hiy = am(? 
for the scattered photons and hv = aimc? for the target photons. We solve eq. 
(HI numerically using the full Klein-Nishina cross section for a single electron [see 



e.g. Illlj. Oimax accounts for the kinematic restrictions on IC scattering. Electrons 
escaping the blob are considered via iesc,rad = ?7^rad/c with rj = 10, which is the 
responding timescale of the electron system. Hence the variability time scale in the 



observer's frame is [see e.g. |l2|: 



, ^csc.rad /r\ 

ivar OC (5) 



To determine the time-dependent model SED of blazars the PDE for the differential 
photon number density, obtained from the radiative transfer equation, is solved 
time-dependently 

^ ^ Rs- ca„Np^[iy, t) + Re f . (6) 

The synchrotron production rate Rs is calculated using the Melrose approximation, 
the inverse Compton rate R^ is treated with the Klein-Nishina cross section, see 
[2| ■ The self absorption of the synchrotron photons is described by Ui^ [2, ll3j . The 
photon loss rate is chosen to be the light crossing time. 
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3. Results 



PKS2155 is used to compare the model to observations. The recently published 
Fermi LAT data [5[ gives a strong constrain on the curvature of the VHE peak in 
the SED of PKS2155. In context of a SSC model this leads to a deep dip between the 
synchrotron and inverse Compton peak in the spectrum. This can only be modelled 
selfconsistently with a more complex electron distribution than a powerlaw used 
in many SSC models. We can produce an electron distribution in the acceleration 
zone with a rising part only by shifting the injection to moderate energies, i.e. 
7o ~ 3 • 10^, which finally is able to explain the observed SED. Table 1 summarizes 
the parameters of the model used to fit the lowstate of PKS2155 shown in Fig. la. 
The model SED is weh in equipartition with h = B^ / {Stt J ^Ned-ymc^) = 0.13. We 



Table 1. Model parameters used to fit the lowstate SED of 
PKS2155. 



Qo(cm-3) 


B{G) 


-Race (cm) 


-Ri-ad(cm) 


f^acc/*csc 


a 


<5 


8.0-105 


1.4 


1.0-1013 


5.0-10" 


1.13 


20 


49 



neglegted the optical data assuming a contribution besides the model region (i.e. 
disk radiation and/or another part of the jet). However the flux level, polarization 
and variability indicate synchrotron radiation from the jet which is considered in 
other models 5]. This might be seen as a major drawback of our fit. This lowstate is 
now used to model the famous VHE outburst of PKS2155 in 2006 [6|. In our model 
this is done by injecting more monoenergetic electrons at 70 into the acceleration 
zone for certain amounts of time having a reasonable physical interpretation as 
fluctuations in the electron density along the jet axis. These electrons also undergo 
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Fig. 1. a) The lowstate of PKS2155 as observed in [5||, red points indicate simultaneous data, 
blue ones are archival |14|. The solid line represents the model SED. The ATOM datapoints have 
been neglegted in the modelling (see text), b) The time averaged SED over the whole flare of our 
model and H.E.S.S. 6] (thick solid line) and the time evolution of the model SED during the flare. 
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acceleration and enter the radiation zone giving rise to the hghtcurve for measured 
flux above 200 GeV in Fig. 2a. The variabihty time scale from Fig. 2a is roughly 
ivar ~ 500 s when concerning the FWHM of the narrowest peak, this corresponds 
to the variability timescale given in eq. ([5]) for the parameters in Table 1. Fig. 
2b shows the time evolution of the electron distribution in the acceleration zone 
during the first peak in the hghtcurve of PKS2155 out of the lowstate with its 
rising and falling powerlaw components and the exponential cutoff. [6| also shows 
the SED of PKS2155 in the high state during the flare as a time averaged spectrum 
over the flare. The result of time averaging the produced model SEDs during the 
flare is shown in Fig. lb. With the model however the time evolution of the SED 
during an outburst is also accessible. Fig. lb. This evolution can not be observed 
directly due to the low flux levels from blazars in the air Cherenkov telescopes, so 
it is important to i) constrain the model with the available data in order to ii) get 
further information about the variability of blazars over the time evolution of the 
SED in the model. 

4. Conclusions 

The presented model is able to motivate the injected powerlaw electron spectra 
often used in onezone SSC models using acceleration and radiative cooling within 
the blob. The spectral electron index arises from the efficiency of the acceleration 
processes compared to the catastrophic losses in the acceleration zone. The posi- 
tion of the high energy cutoff is due to the equilibrium of synchrotron cooling and 
Fermi-I/-II acceleration, while its shape is dominated by the stochastic part of the 
acceleration. By injecting monoenergetic electrons at moderate 70, instead of low 
70S, into the blob our model can explain electron spectra with a rising powerlaw 
part, as suggested in [3|. Therefore the Fermi-II processes are essential. It turns out 
such electron spectra are mandatory to explain the curvature of the VHE peak and 




1000 2000 3000 4000 5000 6000 7000 8000 9000 
t/s (from MJD 53944.0) 




Fig. 2. a) The Hghtcurve of the flux above 200 GeV as measured by H.E.S.S. in 2006 and our 
model, b) Time evolution of the electrons in the acceleration zone during the first peak in the 
hghtcurve shown in a). 
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the deep dip between the two humps in the SED of PKS2155 as constrained by 
the latest Fermi observations. The moderate 70 factor has a physically reasonable, 
but highly speculative, explanation as counterstreaming electrons, since relativistic 
superposition yields 70 = {1-{T^/T'^ - l + "fu\/ll ~ l(r2+72- l)"^)^)"^/^, where 
F is the Lorentzfactor of the blob and 7^ the one of the opposed moving electrons 
entering the blob. Our model is able to fit the multiwavelength data, apart from 
the optical data, in [5| within equipartition (h = 0.13) by shifting the injected elec- 
trons to 7o — 3100. Using this lowstate of PKS2155 as a starting point we are able 
to model the lightcurve of the VHE outburst measured by H.E.S.S. in 2006, even 
without additional parameters compared to simple onezone SSC models. Therefore 
we just have to inject more electrons at 70 into the blob for certain amounts of time 
thus no additional sets of parameters have to be assumed. The time average of the 
model lightcurve fits the H.E.S.S. time averaged SED of the flare. All this leads to a 
selfconsistent and highly constrained description of a HBL and the flaring process. 
With the model we are able to plot time resolved SEDs during an outburst which 
are inaccessible in the VHE regime by direct observations. According to Fig lb the 
changes in the X-ray regime during the flare of PKS2155 are beginning well before 
the H.E.S.S. observation. An average over this time would e.g. lead to a shift of 
the synchrotron peak or a different ratio between the two peaks. Thus one has to 
be carefuU when comparing SEDs especially when they have been averaged over 
different times. To make these predictions quantitatively more exact, simultaneous 
observational data during an outburst is required. 
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